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ABSTRACT 

Ram pressure stripping can remove significant amounts of gas from galaxies that 
orbit in clusters and massive groups, and thus has a large impact on the evolution 
of cluster galaxies. In this paper, we reconstruct the present-day distribution of ram- 
pressure, and the ram pressure histories of cluster galaxies. To this aim, we combine 
the Millennium Simulation and an associated semi-analytic model of galaxy evolution 
with analytic models for the gas distribution in clusters. We find that about one 
quarter of galaxies in massive clusters are subject to strong ram-pressures that are 
likely to cause an expedient loss of all gas. Strong ram-pressures occur predominantly 
in the inner core of the cluster, where both the gas density and the galaxy velocity are 
higher. Since their accretion onto a massive system, more than 64 per cent of galaxies 
that reside in a cluster today have experienced strong ram-pressures of > 10~ n dyn 
cm~ 2 which most likely led to a substantial loss of the gas. 
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1 INTRODUCTION 

In clusters, galaxies can lose some or all of their gas by ram 
pressure stripping (RPS) due to their motion through the 
intracluster medium (ICM). Both analytical estimates and 
hydrodynamical simulations show that RPS can remove a 
significant amount of gas from galaxies, and can thus explain 
observations such as the HI deficiency of cluster disc galax- 
ies (see e.g. lHavnes fc Giovanellil I l986l ; ISolanes et al.ll200ll; 
ICavatte et all 19941 ), and the truncated star forming discs 



in th e Virgo cluster ijWarmelsl 1 19881 : iKoopmann fc Kennevl 
120041 ). 

RPS only affects the gaseous component of the galaxy 
so that a distinct feature of ram pressure stripped galaxies 
is that their gas discs are distorted or truncated while their 
stellar discs appear undisturbed. An increasing number of 



available in the last years (see for example Irwin et al. 


1987 


Veilleux et al.lll99S 


, Kennev et al.ll2004,IVollmer et al. 


2004 


Chung et al.ll2007l. 


Sakelliou et al.ll2005l; ISun et alj|2007l). 



RPS is commonly cited in early work as a possible 
explanation for the increased fraction of SO galaxies in 
rich clusters relative to t he field l|Biermann fc Tinslevl[l~975l 
iButcher fc Oemlerl This explan ation was dismissed 

in the original paper bv lPressled (|l980l ) on the basis of the 
observation that the relation between different galaxy pop- 
ulations and local density appears to hold independently of 
the cluster richness. Later studies pointed out that addi- 
tional mechanisms that lead to a significant redistribution 
of mass and/or new star formation are required to explain 



the entire SO pop ulation of galaxy clusters (see for example 
iMoran et al.ll2007l and references therein). 

The role of RPS in the chemical enrichment of the 
ICM has also been discussed. Observational data suggest 
that the ICM is enriched with metals to approximately one 
third of the solar value, suggesting that some fraction of 
the metals must have originated from cluster galaxies and 
since been removed from them. Processes responsible for the 
supply of this enrich ed gas include AGN feedback (see e.g. 
iRoediger et alJfeoOTl and references therein), galactic winds 
driven by supernovae explosions, and ram-pressure str i pping 
ll White! fl99ll; iMori fc Burkerd [2OO0I : ISchindler et~ai] (20051 ; 
iDomainko et al.ll2006l ). It should be noted, however, that 
although RPS certainly affects the metallicity of the ICM, 
it may n ot be the dominant me chanism. Numerical simula- 
tions by IDomainko et al.l l|2006f ) indicate that RPS can ac- 
count for only about 10 per cent of the ICM metal content 
within a radius of 1.3 Mpc. 

The first analytical estim ate of RPS dates back to the 
paper bv lGunn fc Gottl (Il972l ) who proposed that for galax- 
ies moving face-on through the ICM the success or failure of 
RPS can be predicted by comparing the ram pressure with 
the galactic gravitational restoring force per unit area . Later 
hydrodynamical s i mulations of RPS dAbadi et al. 19991. 
Q/uilis et ail |2000|. ISchutz fc Struck l200ll. iMarcolini etall 
2003!; lAcreman et all 120031, IRoediger fc, Briiggen! 120061 . 
Roediger et all 120061 . IRoediger fc Briiggen! 120071 ') suggest 



that this analytical estimate fares fairly well as long as the 
galaxies are not moving close to edge-on. The ICM-ISM in- 
teraction is, however, a complex process influenced by many 
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parame ters. Different aspect s have been studied by several 
groups. iRoediger fc Henslerl (|2005h and others have shown 
that the ICM-1SM interaction is a multistage process: The 
most important phases are the instantaneous stripping, on 
a time-scale of a few 10 Myr, an intermediate phase, on a 
time-scale of up to a few 100 Myr, and a continuous strip- 
ping phase that, in principle, could continue until all gas is 
lost from the galaxy. 

While numerical simulations and observations indi- 
cate that RPS has important consequences on the amount 
of gas in cluster galaxies, this physical process is usu- 
ally not included in semi-analytic models of galaxy for- 
mation. The effect of ram-pressure stripping has been dis- 
cussed only in a couple of studies using semi-analyti c 
models (jOkamoto fc Nagashimall2003l : lLanzoni et al.ll2005h . 
and is shown to produce only mild variations in galaxy 
colours and star formation rates. This happens because 
the stripping of the hot gas from galactic haloes (stran- 
gulation) suppresses the star formation so efficiently that 
the effect of ram-pressure is only marginal. We note that 
the studies mentioned above include ram-pressure strip- 
ping based on the analytical criterion formulated origi- 
nally in Gunn fc Gott (19 72). In recent numerical work, 
IRoediger fc Briiggenl i|2007l ) have shown that this formula- 
tion often yi elds incorrect mass lo ss rates. Other numerical 
studies (e.g. IVollmer et al.l l^OPl! )) have argued that ram- 
pressure stripping can also temporarily enhance star forma- 
tion. An updated modelling of ram-pressure stripping that 
takes into account these results has not been included yet 
in semi-analytic models. We plan to address this in future 
studies. 

For a study of the ram pressure distribution, it is nec- 
essary to have information on the dynamics of galaxies and 
on the properties of the ICM. The dynamics of dark mat- 
ter halos has been studi ed in a number of papers using nu- 
merical simulations (e.g. Benson 2005; Khoc hfar fc Burkertl 
120051 ; [Diemand et al.ll2004T ). However, we know of no study 
on the distribution and history of ram pressures experienced 
by galaxies in clusters. If ram-pressure plays some role in es- 
tablishing the observed morphological mix in galaxy clusters 
and/or the observed radial trends, it is important to quantify 
the distribution and history of ram pressures experienced by 
galaxies that reside in clusters today. This is the subject of 
this paper. 



2 METHOD 

For this s tudy, we rely o n the Millenium simulation de- 
scribed in ISpringel et al.l l|2005l ). This largest dark matter 
simulation to-date follows N — 2160 3 particles of mass 
8.5 x 10 8 h~ l M & within a comoving box of size 500 h~ x Mpc 
a side. The underlying cosmological model is a ACDM model 
with fl m = 0.25, fl b = 0.045, Q A = 0.75, n = 1, (T 8 = 0.9 
and h = 0.73, where the Hubble constant is parametrised 
as Ho — 100 h km s _1 Mpc _1 . Given its high resolution 
and large volume, this simulation allows us to make statisti- 
cally significant inferences about the ram pressure histories 
of galaxies in a representative sample of clusters. 

During the simulation, 64 snapshots were saved, to- 
gether with group catalo gues and their e mbed ded substruc- 
tures. As explained in ISpringel et al.l (|2005h , dark mat- 



ter haloes are identified using a standard friends-of-friends 
(FOF) algorithm with a linking length of 0.2 in units of 
the mean particle separation. Each FOF group is then de- 
composed into a set of disjoint substructures identified as 
a locally overdense region in the density field of the back- 
ground halo. The selfbound part of the FOF group itself 
will then also appear in the substructure list and represents 
what we will refer to below as the main halo. This particu- 
lar subhalo typically contains 90 per cent of the mass of the 
FOF group. The group catalogues were then used to con- 
struct detailed merging history trees of all gravitationally 
self-bound dark matter structures. These merger trees form 
the basic input needed by the semi-analytic model used in 
|Pe Lucia fc BlaizotJ (|2007f ). 

We extracted the orbital parameters of the galaxies 
from the public archive of the Millenium Run data bascQ, 
and we refer to the original paper for details about the phys- 
ical processes that are part of the model. Our analysis uses 
only the orbital parameters of model galaxies and therefore 
does not rely explicitely on the details of the semi-analytic 
model itself. One important limitation to take into account 
is that most of the galaxies in massive clusters are " orphan" 
galaxies, i.e. galaxies that are not as sociated with dark 
matte r substructures. As explained in |Pe Lucia fc BlaizotJ 
(2007) and in previous papers related to this model, sub- 
structures allow us to trace the motion of the galaxies sit- 
ting at their centres only until tidal truncation and strip- 
ping disrupt the substructures down to the resolution limit 
of the simulation (whi ch for the Millennium Simulation is 
L7x 10 10 M h _1 ) (e.g. lDe Lucia et alJliooi iKravtsov et al.l 
|2004| ). After this time, the galaxy is assumed to merge onto 
the central galaxy of its own halo on a dynamical friction 
time-scale, and its position and velocity are traced by track- 
ing the position and velocity of the most bound particle 
of the halo at the last time there was a substructure. As- 
suming that the position and velocity of the most bound 
particle at the last time the substructure could be identified 
serve as correct initial conditions to track the orbit of the 
galaxy, t he ensuing ram pres sure will also be correct. Recent 
studies (jConrov et alj 120071 ) have argued that a significant 
fraction of the satellite population from disrupted subhaloes 
is unbound and goes to the intra-cluster light component. A 
large fraction of galaxies in massi ve haloes is represent ed by 
orphan galaxies. If, as argued in IConrov et al.1 l|2007l ). the 
model we have used in our study leaves behind an excess of 
orphan galaxies, this would affect some of the results pre- 
sented in this study. However, the issu e regarding o rphan 
galaxies does not seem to be settled. IWang et al.l (|2006t ) 
have shown that orphan galaxies are needed in order to re- 
produce the observed correlation function on small scales. 
The existence of intra-cluster light suggests that tidal ef- 
fects or mergers can unbind some of the stars in the satellite 
galaxies. Published results, however, offer little indication 
of appropriate recipes for treating this process within semi- 
analytic models. Observationally, the total amount of the 
intra-cluster light is very difficult to estimate and published 



1 A description of the publicly available catalogues, and a 
link to the database can be found at the following webpage: 
http://mpa-garching.mpg.de/millennium/ 
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esti mates vary from less than 20% to more than 50% (see 
e.g. IZibetti et af1l2005l ; iGonzalez et"alll2005h . 

As the Millenium run is a dark-matter only simulation, 
we have to make assumptions about the distribution of the 
gas in order to compute the ram pressure exerted on the clus- 
ter galaxies. A first estimate is to app roximate the ICM as 
isothermal and hydrostatic in a NFW (jNavarro et al.|[l996l ) 
halo whose profile is given by: 



pDM(r) 



ScPcO 



(1) 



(r/r s )(l + r/r s ) 2 ' 
where p c o is the critical density of the universe at z = 0, and 



S C (M) « 
r s {M) = 



3 x 10 3 n [l + z t (M)} 3 , 



1 / 3M 



ry ir (M) _ 

c(M) ~~ c(M) \4irA c p c0 



1/3 



(2) 
(3) 



In the above expressions, £Iq is the density parameter at z = 
0, A c (r2o,Ao) is the collapse factor in a spherical nonlinear 
model, Zf(M) is the average formation redshift of objects of 
mass M, and c(M) the concentration parameter. 

For an isothermal spherical gas cloud with temperature 
Tx, the density distribution p g in hydrostatic equilibrium 
satisfies the equation: 



fcTx d In p g 



GM(r) 



(4) 



/im p dr r 2 

where p and m p denote the mean molecular weight (we 
adopt 0.59 below) and the proton mass. If one neglects the 
gas and galaxy contributions to the gravitational mass in 
the right-hand side, then the mass enclosed within a radius 
r, can be obtained from Eq. {1} and is given by: 



M(r) = 4n5 c p c0 r 3 .m{r/r s ), 



(•») 



where m(r/r s ) is the function m evaluated at r/r s and m is 
given by 



m(x) = ln(l + x) 



1 + x 



Equation Q can be integrated analytically to give 

ln(l + r/r s ) s 



ft! 00 



with 



Pgo exp 



27b 



(6) 

(7) 
(8) 

27fcTx ' (9) 

as shown in iMakino et al.l (|l998l ). The cluster gas temper- 
ature Tx is expected to be close to the virial temperature 
T v i r (M) of the dark matter halo. In the profile {TJ, the latter 
is in fact dependent on the radius r: 



2 V' r/r s 
p s0 exp(-27fe/2) (l + r/r s ) 27b/(2 ^ ) , 



8iYG(j.m p 6 c (M)p c0 r 2 



kT vil (r) = 7 



Gpm p M (r) 
37 : 



(10) 



where 7 is a fudge factor of order unity which should be 
determined by the efficiency of the shock heating of the gas. 



lEke" et al.l l|l998l ) adopted 7 = 1.5 as their canonical value 
in the analysis of X-ray cluster number counts, and this is 
also the value we use here. Substituting (|10p into equation 
([9|, one finds that: 



Mr) = 2lL. 

y ' 97 r s 





r 








r + r s _ 



(11) 



For an absolute determination of the gas profile, we 
also need an estimate of the cluster gas fraction which 
we chose to be equal to the universal gas mass fraction 
/gas = tt b /tt m = 0.022ft- 2 /0.3 = 0.14. There is evidence 
that in massive structures such as gal axy clusters this 
gas fraction is fairly constant over time (| Allen et al.ll2007l : 
iLaRoque et al.l 2006). The ram-pressures computed in the 
following scale simply with this gas fraction. 



As shown by iNavarro et al.l l|l996l ) , at large radii the 
density profile of an isothermal gas drops less rapidly 
than the dark matter (see their Fig. 14). This levelling 
(which is observed out side a radius ~ 2 x -R v j r ) is not ob- 
ser ved in real clu s ters (Finoguenov et al.ll200ll ). The model 
by IMakino et al.l l|l99£ ) has been later extended to non- 
isothermal gas with a polytropic eq uation of state (e.g. 
ISuto et al.lli998l ; Ascasibar et al. 2003, IVoit et~ai]|2002l ). 

One alternative approach t o reconstruct clust e r mas s 
profiles has been suggested by iKomatsu fc Seliakl (|200 lh - 
They present an analytic method to predict gas density and 
temperature profiles in dark matter haloes that does not 
rely on the isothermal approximation. In this model, the 
gas density profile traces the dark matter density profile in 
the outer parts of the haloes (an assumption that is also sup- 
ported by hydrodynamic simulations), and the gas obeys a 
polytropic equation of state. In the inner regions of galaxy 
clusters, gas temperature often increases with radius up to 
100 — 200 kpc and then mildly decreases in the outer regions. 
The additional assumption that the gas temperature has to 
vary monotonically with density therefore limits this model 
to regions outside th e inne r 100 — 200 kpc. In the model 
bv lKomatsu fc Seliakl (|200ll ). the gas density distribution is 
given by: 



pgas(r) = p g a S (0)y g as(r/r s ), 
where p gas (0) is the gas density at r — and 



(12) 



7-1 



(x) 



m(c) 



[l-ln(l+a;)/a;] , (13) 



where again we assume a NFW density profile of concentra- 
tion c. The effective value for 7 is found to depend weakly 
on the concentration parameter according to the following 
equation: 



7 = 1.15 + 0.01 (c- 6.5), 
and the parameter r\ is given by 



(14) 



77(0) = 0.00676 (c - 6.5) 2 + 0.206 (c - 6.5) + 2.48. (15) 
Equations (|12[) - (|15[ ) allow a reconstruction of the gas profile, 
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and, with the velocity information of the galaxy, of the ram 
pressure. 

Both models discuss ed above, the isothermal and the 
iKomatsu fc Seliakl l|200lf) model, rely on the assumption 
of hydrostatic equ ilibrium. Numerical simulations (e.g. 
lAscasibar et al.ll2003T ) show that that this is not a bad ap- 
proximation unless they have suffered a major merger in „ 
their recent past. However, there is some indication that the g 
dynamic ICM can lead to variations in the ram pressure. g. 
E.g. it has been suggested that in the Virgo cluster sloshing c£ 
motions of the intracluster gas lead to changes in the ram : 
pressure l|van Gorkomll2007r i. \ 



3 RESULTS 

3.1 Present-day distribution of galactic ram 
pressures 

In this section, we discuss the distribution of ram-pressures 
of cluster galaxies at the present epoch (i.e. at z = 0). For 
this analysis, we have selected a number of massive haloes 
from the Millennium data base, and identified all the galax- 
ies from the De Lucia & Blaizot (2007) semi-analytic cata- 
logue that share the same FOF group. For this analysis we 
have excluded the central galaxy of each FOF group (i.e. the 
brightest cluster galaxies). 

Our sample is composed of 20 clusters with masses 
close to M = 10 15 M© (M 20 o between 7.3 • 1O 14 M and 

1.2 ■ 1O 15 M ) and 174 clusters with masses close to M = 
10 14 M Q (M 200 between 9.7 ■ 1O 13 M and 1.03 ■ 1O 14 M ), 
yielding a total of 78,178 and 74,294 galaxies, respectively. 

Fig. [T]shows the distribution of instantaneous ram pres- 
sures exerted on galaxies within clusters of masses M — 
1O 14 M (black histogram) and M = 1O 15 M (red his- 
togram). For the M = 1O 14 M clusters, the distribution 
peaks at ~ 10" 13 dyn cm" 2 , while for the M = 1O 15 M 
clusters, it peaks at ~ 10 -12 dyn cm -2 (see caption of Fig.[T] 
for exact values of the median and the mean). The shapes 
of the distributions for both mass ranges are very similar, 
with a tail at lower ram pressures and a fairly sharp cutoff at 
hi gher ram pressures. The co rresponding plot for the model 
by IKomatsu fc Seliakl (|200ll ) is shown in Fig. [2] The peaks 
of the two distr ibutions are at sim i lar ra m pressure values. 
However, in the lKomatsu fc S eliak ( 200l|) model the cut-off 
at high ram pressures is less sharp, and the distributions 
appear to be less skewed. This is because in the isothermal 
model, the density is underestimated at large radii. Thus 
galaxies at large distances from the centre (which make up 
the tail of low ram-pressure values that is visible in Fig. 1) 
suffer a lower ram pressure in the isothermal model com- 
pared to the Komatsu model. In both Figs. 1 and 2, the 
solid histograms show the mean obtained for all the clusters 
in each mass bin considered, while the error bars indicate 
the scatter of the distributions. 

According to GG72, the mass lost from a galaxy 
through RPS depends on the gravitational restoring force 
per unit area. Clearly, other factors, such as inclination, gas 
content, morphological type, star formation rates, etc. also 
play a potentially significan t, albeit poorly unde r stood , role. 
The catalogue compiled by |Pe Lucia fc Blaizod (2007) con- 
tains information that would allow us to estimate the mass 
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Figure 1. Distribution of galactic ram pressures in isothermal 
clusters with virial mass M = 1O 14 M (black) and M = 1O 15 M 
(red) at z = 0. For M = 1O 15 M , the mean ram pressure is 



10 



-10.8 



dyn cm 2 and the median 10 118 dyn cm 2 . For M 



10 Mq, the mean ram pressure is 10 dyn cm -2 and the 

median 10 -12 ' 7 dyn cm -2 . The error bars denote the cluster-to- 
clustcr scatter. 



30 - 



°g(p, m) 

Figure 2. Distribution of galactic ram pressures in clusters with 
mass M = 1O 14 M (black) and M = 1O 15 M (red) at z = 
for the model by Komatsu & Seljak (2001). For M = 1O 15 M , 
the mean ram pressure is 10 — 10,7 dyn cm -2 and the median ram 
pressure is 10 — 11-9 dyn cm -2 . For M = lO 14 Af , the mean ram 



pressure is 10 



dyn cm and the median ram pressure is 



dyn cm 
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loss due to RPS. However, ram-pressure stripping is not self- 
consistently included in the model. Therefore, we focus here 
on merely computing the ram-pressures experienced by the 
galaxies at a given time a nd position within the clust er. 

In the simulations bv lRoediger fc Briiggenl l|2006t ). ram 
pressures of ~ 1CP 12 dyn cm -2 were called weak, ~ 1CP 11 
dyn cm -2 medium and ~ 10 -10 dyn cm -2 strong. Subject 
to strong ram pressure, a typ ical spiral galaxy as simulated 
in lRoediger fc Briiggenl (120061 ) with mass ~ 2 ■ 10 n M Q loses 
all its gas within ~ 50 Myrs. Medium ram pressure removes 
approximately half of the gas within ~ 200 Myrs, and weak 
ram pressure removes relatively little gas. These numbers 
depend of course on the structure of the gaseous, stellar and 
dark matter component, and should just serve for orienta- 
tion. Figs. [1] and [2] show that in a M = 10 15 M Q cluster, 
approximately 27 per cent of galaxies experience ram pres- 
sures of > IP" 11 dyn cm 2 in the isothermal ICM model. 
In the iKomatsu fc Seliakl (|200lh model, the corresponding 
fraction is about 24 per cent. In a M = 10 14 Mq cluster, 
these numbers are 10 % for the isothermal model and 9 % 
for the IKomatsu fc Seliakl (|200lh model, respectively. 

Since ram-pressure depends on the density of the gas 
and on the velocity of the galaxy, it is expected to be 
stronger closer to the centre. In Fig. [3] we plot the ram pres- 
sure as a function of radius for the Komatsu (2001) model 
and a 10 15 M Q mass cluster. The distribution shows a sharp 
upper edge which is determined by the escape velocity at this 
radius and the gas density at that position. In the isothermal 
model, the upper envelope can be approximated by: 



max / \ _ 
Pram M » /O, 



2GM(r) 



(16) 



If we approximate M(r) by Eq. ([5]) and p g by Eq. (|12p . we 
can rearrange the result to give 



-2B(x) 

pSmOc) « 87vGp sQ 5 c p c0 rs — r (1 + x) 2B{x)/x , (17) 



B(x) 



where x — r/r B and 



B(x) = x[ln(l +x)-x/(l + x)\ 



(18) 



Eq. (|17[1 describes the maximum ram pressure pretty 
well for x > 0.2. Similar, though less simple, expres- 
sions can be found for t he maximum ram pressure in the 
IKomatsu fc Seliakl (l200ll) model. For the two gas models, 
the isothermal and the Komatsu fc Seliakl (|200ll ) model, the 
distributions are somewhat different. While within the virial 
radius the ram pressures are very similar, at larger radii, 
the isothermal mo d el yie lds lower ram pressures than the 
IKomatsu fc Seliakl (l200ll ) model. The latter model allows 
the temperature to decrease at larger radii which leads to a 
higher density than in the isothermal model. This also ex- 
plains the narrower distribution of ram pressures shown in 
Fig. [2]with fewer galaxies in the low ram pressure tail of this 
distribution. 



3.2 Ram pressures histories 

In this section, we analyse the ram-pressure history suffered 
by galaxies that reside in a cluster today. In the following, 



0.4 0.6 

distance from centre/ rvir 



Figure 3. Distribution of galactic ram pressures (in dyn cm 2 ) 
with radius in a cluster with mass M = 1O 15 M at z = (model 
by Komatsu & Seljak 2001). The thick solid line denotes the mean 
of the distribution. The dashed line shows the mean of the distri- 
bution for the isothermal model. The yellow points indicate the 
ram pressures of the individual galaxies in the Komatsu & Seljak 
model. 



we restrict our analysis to galaxies with with a B-band mag- 
nitude mags < —19, where mags is given in the De Lucia 
& Blaizot (2007) catalogue. For each galaxy, we walk its 
merger tree backwards in time by linking it with its most 
massive progenitor at each snapshot until the galaxy be- 
comes a central galaxy of a FOF halo. In Fig. [4] we plot 
a random selection of the ram pressure histories of galaxies 
that end up in M = 10 15 Mq clusters. This and the following 
plots re fer to an isotherma l ICM (the corresponding results 
for the IKomatsu fc Seliakl l|200ll) model are very similar). 
Fig. 2] shows that the ram pressures fluctuate strongly with 
time. No strong trend with redshift is visible, showing that 
galaxies underwent phases of strong ram pressure even at 
high redshifts. In Fig. [5] we show how the ram pressure (top 
left panel), the mass of the parent FOF group (top right 
panel), the relative velocity and the distance to the cluster 
centre (bottom left panel), and the ambient ICM density 
(bottom right panel) vary as a function of redshift for a ran- 
domly selected galaxy that reside in a M = 10 15 M Q clus- 
ter at the present epoch. This plot shows that the galaxy 
attains high velocities when it passes close to the cluster 
core. In this region, the ambient ICM density is also high- 
est, such that the highest ram pressures values are obtained. 
We note that the outputs of the simulation are not sampled 
finely enough in time to reconstruct the galaxy ram pressure 
history very accurately. The vertical lines in the top right 
panel of Fig. [5] show the redshifts of the simulation output. 
The time between snapshots is too large to sample the or- 
bits of the galaxies with great precision. In some cases, the 
ram pressure fluctuates by 1-2 orders of magnitudes between 
snapshots and it is difficult to assess in each case what hap- 
pens in between. Consequently, the ram pressure histories 
give conservative bounds on the maximum and minimum 
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Figure 4. Ram pressure histories of a random sample of galaxies 
that end up in a M = 1O 15 M0 cluster. 



ram pressures suffered by each galaxy in the course of its 
life. We also note that periods of strong ram pressure are 
often interspersed with periods of weak ram pressure. 

In these periods, the interstellar medium could be ac- 
creted back onto the galaxy or replenished (by e.g. merg- 
ers with gas-rich satellites that have not yet suffered sig- 
nificant stripping). The median redshift at which a galaxy 
has most recently experienced moderate or strong ram pres- 
sure is close to 0.1 for all galaxies in our sample. Strong 
ram-pressure episodes are, however, expected to have a sig- 
nificant effect on the following evolution of the galaxy. In 
Fig. [6] we show the maximum ram-pressure that a galaxy 
(or its most massive progenitor) has experienced since the 
time of accretion. The solid histogram shows results for a 
M — 1O 14 M0 cluster, while the dashed histogram refers to 
a M = 1O 15 M cluster. The fi gure shows that in a massive 
cluster, more than 64 per cent of galaxies have experienced 
ram pressures of > 10 -11 dyn cm" 2 , and 32 per cent of 
galaxies have had ram pressures greater than > 1CP 10 dyn 
cm~ 2 . Ram pressures of this magnitude most likely strip 
the galaxy of all its gas in a short time interval of a few mil- 
lion years. The corresponding fraction for a M = 1O 14 M0 
cluster are lower but still significant (52 per cent and 11 
per cent respectively). If ram-pressure stripping is responsi- 
ble for the morphological transformation of spiral galaxies 
infalling onto the cluster from the field, these numbers in- 
dicate that about half of the galaxies residing today in a 
M = 1O 14 M0 cluster, and a larger fraction for more mas- 
sive clusters, should be gas poor. Ellipticals and lenticulars 
make up about 70 — 80 per cent of the galaxy population 
of massive clusters in the local Universe (see e.g. Dressier 
1980). A significant fraction of these could be therefore en- 
tirely explained by ram-pressure stripping. 

The original work by Dressier (1980) also showed that 
there is a clear trend for an increasing fraction of early 
type galaxies with decreasing distance from the cluster cen- 
tre. Fig. [7] shows the fraction of galaxies that have suffered 
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Figure 5. Properties of one selected galaxy versus redshift. Top 
left: Ram pressure (in dyn cm -2 ). Top right: Mass of the most 
massive FOF group (i.e. the host cluster or group) of this galaxy. 
Vertical lines mark the redshifts of the simulation output. Bottom 
left: Velocity (black, solid line) and distance to cluster centre (red, 
dashed). Botton right: Ambient ICM density (using the Komatsu 
k, Seljak (2001) model. 
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Figure 6. Maximum ram pressure that a galaxy from a M = 
10 14 M Q (solid line) and M = 10 15 M Q cluster (dashed line) or 
its most massive progenitor has experienced during its life time. 
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strong (> 10 -10 
and weak (> 10 
of their current (z = 0) distance from the cluster centre. 
The figure clearly shows that a larger fraction of the galax- 
ies that reside in the cluster core have suffered significant 
ram-pressure. This fraction monotonically decreases with 
distance from the cluster centre, in qualitative agreement 
with the observed trends. In M — 10 15 M Q clusters, virtu- 
ally all galaxies in the inner 300 kpc have suffered strong 
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restricted to within the virial radius. Our results show, 
however, that virtually all the galaxies in clusters suffered 
weaker episodes of ram pressure, suggesting that this phys- 
ical process might have a significant role in shaping the ob- 
served properties of the entire cluster galaxy population. 

The limited number of simulation outputs does not 
allow us to reconstruct accurately the orbit of the clus- 
ter galaxies, and therefore their ram-pressure histories. 
Our result show that ram pressure fluctuates strongly so 
that episodes of strong ram-pressure alternate to episode 
of weaker ram-pressure. During these time intervals, the 
gaseous reservoir could be replenished and new episodes of 
star formation could occur. Our results indicate that ram- 
pressure stripping must play a significant role in the evo- 
lution of galaxies residing in massive clusters. A more self- 
consistent modelling is however required in order to draw 
more quantitative conclusions about the importance and ef- 
fects of this physical process. 



dist. from centre (Mpc) 



Figure 7. Fraction of galaxies that have suffered strong (solid 
line), medium (long-dashed line) and weak (short-dashed line) 
ram pressure in the course of their history as a function of their 
current (z = 0) position in the cluster. The black lines correspond 
to M = 10 15 M Q clusters and the red lines to M = 10 14 M ( 
clusters. 







ram pressures after accretion. For M = 1O 14 M0 clusters, 
this fraction decreases to about 2/3. Nearly all galaxies in 
clusters have experienced medium ram pressures and have 
thus been influenced to some degree by ram pressure strip- 
ping. 
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4 CONCLUSION 

In this study we have taken adva ntage of the Millennium 
simulation bv lSprineel et all (2005) and of the publicly avail- 
able semi-analytic model by De Lucia & Blaizot (2007) to 
reconstruct the ram-pressure distribution and ram-pressure 
history of galaxies that reside in clusters at the present 
epoch. The gas profile in dark matter is described through 
two analytic models which give fairly similar results. We 
have not included ram-pressure stripping self-consistently 
in the semi-analytic model employed for our study. Instead 
we have simply used the available information about the or- 
bital distribution and galaxy merging trees to estimate the 
importance of ram-pressure stripping on galaxies that reside 
in massive haloes at the present epochs. 

We find that more than half of the galaxies in a M = 
1O 15 M0 cluster, have experienced ram pressures of > 10 -11 
dyn cm -2 after their accretion onto a massive system. This 
fraction is only slightly lower for M = 10 14 Mq clusters, im- 
plying that a significant fraction of galaxies in clusters at the 
present epoch suffered substantial gas loss due to ram pres- 
sure stripping. The fraction of galaxies that suffered signif- 
icant ram-pressure after accretion increases with decreasing 
distance from the cluster centre, in qualitative agreement 
with the observed increase of early- type galaxies. 

As expected, strong episodes of ram-pressure occur pre- 
dominantly in the inner core of galaxy clusters, and are 
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